One-dimensional Dirac oscillator in a thermal bath 
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We analyze the one-dimensional Dirac oscillator in a thermal bath. We found that the heat 
capacity is two times greater than the heat capacity of the one-dimensional harmonic oscillator for 
higher temperatures. 
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The name Dirac oscillator was first introduced by Moshinky and Szczepaniak |l| for a Dirac equation in which 
the momentum p is replaced by p — > p— imuftr with r being the position vector, m the mass of the particle and u) 
the frequency of the oscillator. The name was given because the determination of the spectrum and eigenstate of a 
Dirac oscillator [lj was obtained from that of an ordinary oscillator with a very strong spin-orbit coupling term. The 
concept, though not the name, had been discussed previously first by Ito et al. and later by Cook Q and Ui and 
Takeda and possibly others. Besides intrinsic mathematical interest the study of the Dirac oscillator has invoked 
much attention because of its various physical applications. For instance, the last three chapter of a book of Moshinsky 
and Smirnov p| deal with the Dirac oscillator relativistic interactions between systems of one, two and three particles. 
Moreno and Zentella showed p that it could be the object of an exact Foldy-Wouthuysen transformation. Bcnitez et 
al. found the electromagnetic potential associated with the Dirac oscillator, and showed that this exactly soluble 
problem has a hidden supersymmetry, which is responsible for the special properties of its energy spectrum. They 
also calculated the related superpotential and discussed the implications of this supersymmetry on the stability of the 
Dirac sea. 

Rozmej and Arvieu § have shown a very interesting analogy between the relativistic Dirac oscillator and the Jaynes- 
Cummings model. They showed that the strong spin-orbit coupling of the Dirac oscillator produces the entanglement 
of the spin with the orbital motion similar to what is observed in the model of quantum optics. 

More recently Nogami and Toyama and Toyama et al ^3] have studied the behaviour of wave packets of the 
Dirac oscillator in the Dirac representation in (1+1) dimensions. The aim of these authors was to study wave packets 
which could possibly be coherent. This reduction of the dimension was brought as an attempt to get rid of spin effects 
and to concentrate on the relativistic effects. 

As a matter of fact, the one-dimensional relativistic Dirac equation has been largely used to treat physical problems 
where relativistic effects could play an important role. Indeed, this treatment gives enlightening information about 
more realistic models in solid state physics |l II . particularly in semiconductor theories |l2(. 

In spite of the great number of papers that has been recently published concerning the solution and properties of 
the Dirac oscillator, as far as we know no one has reported on its thermodynamics properties. In order to overcome 
this lack of information, in this letter we study the one-dimensional Dirac oscillator in a thermal bath. 

Let us consider the energy spectrum for a one-dimensional Dirac oscillator |13| 



E n = y/[2(n + l)]huomc 2 + m 2 c 4 , (1) 

where m is the mass of the particle, w is the angular frequency of the oscillator and n is a positive integer number. 
As we can see, for n + 1 << mc 2 /27kJ, the spectrum of the one-dimensional Dirac oscillator is approximated, up to 
Eq, to the spectrum of the one-dimensional harmonic oscillator. 

The partition function of the Dirac oscillator at temperature T is obtained through the Boltzmann factor, 



z= ^ e -(B n -Bb)/9 = e V5/J^ e -V5S+T/J s (2 ) 

71=0 n=0 

were j3 = 1/kT, k the Boltzmann constant, a = 2hujmc 2 , b — m 2 c 4 -I- 2fiiomc 2 , and we are using a frame in the heat 
bath. We have subtracted the ground state energy in order to make a comparison with the non-relativistic harmonic 
oscillator. All thermodynamics quantities for the Dirac oscillator are obtained through the partition function yj. 



* e-mail: renan@fisica.ufc.br 
t e-mail: carlos@fisica.ufc.br 



2 



Let us first analyze the convergence of the partition function. The function f(x) = exp(— \/ax + 6/3) is a monoton- 
ically decreasing function and the corresponding integral 

1(13) = / e -^+ S/3 d a; = — e-^(l + /3Vb), (3) 
Jo «p 

is convergent. Thus from the theorems of convergent series, this implies that the partition function is convergent. 
In order to evaluate the partition function, we can use the Euler-MacLaurin formula 

E /(») = o/(°) + / f to** - E T^TT^/^-^CO), (4) 
n=0 ^ J ° p=l {ZP) - 

where Bi v are the Bernoulli numbers, B^ — 1/6, B4 = —1/30, .... They are defined through the series 

t ^ t n 

B n -. (5) 



e* — 1 ^— i n\ 

n— 

Then the partition function can be written as 

Z = \ + ^(1 + M) - E ^/^(O)' (6) 

To compute the partition function, we need to calculate the sum in the above expression. For our case it can be 
done only by numerical methods. Up to p — 2 this sum can be written as: 

2 B 2 



e ^Eg|/ (M ( ) = 

1 [3 a 1 /3 a 3 1 f3 2 a 3 1 (3 3 a 3 
24 IJt ~ 1920 W 1 ~ 1920 ~~ 6 1 5760 6 3 / 2 



(7) 



Before computing numerically the partition function, let us analyze the case of higher temperatures. This correspond 
to j3 <C 1. All terms in the sum of I© have a positive power in (3, which are very small compared with the term 
2(1 + {3Vb)/a{3 2 . Hence, for (3 < 1 



( g ) 

a(3 z 



Now we can easily obtain the mean energy and the heat capacity of the Dirac oscillator for higher temperatures 



Cc — ^2k. (10) 

These results show that, for higher temperatures, the mean energy and the heat capacity for the Dirac oscillator is 
two times the mean energy and the heat capacity of the non-relativistic one dimensional harmonic oscillator. This is 
due, only to relativistic effects, since in one dimension there is no spin coupling. 

Now we briefly discuss our numerical results on the calculation of the partition function © . First of all, we should 
mention that, using ©, the curves for the three thermal functions, namely, mean energy, heat capacity and free 
energy, are identical to those obtained for the nonrelativistic one-dimensional harmonic oscillator, in the limit of low 
temperatures (0 to A00K). 

We display a comparison between the Dirac oscillator and the harmonic oscillator only, for high temperatures 
(Fig- ffl and Fig. EJ) . It is seen that the free energy and the mean energy are greater for the Dirac oscillator than for 
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FIG. 1: The free energy for the harmonic and the Dirac oscillator for higher temperatures and Tiw/mc 2 = 1. Dashed line 
corresponds to the Dirac oscillator and the solid line the harmonic oscillator. 
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FIG. 2: The mean energy for the harmonic and the Dirac oscillator for higher temperatures and hca/mc 2 = 1. Dashed line 
corresponds to the Dirac oscillator and the solid line the harmonic oscillator. 

the harmonic oscillator. From Fig. [21 it is also seen that the heat capacity for the Dirac oscillator is two times the 
heat capacity for the harmonic oscillator, result that was anticipated by the analytical calculations presented above. 

It is worthwhile to mention that numerical calculations of the partition function for high temperatures (/3 <tc 1), force 
us to handle several very small variables, which could lead to deceiving results. So, it is required a precise estimation 
of the involved physical quantities. Therefore, we choose the angular frequency of the oscillator as oj — 10 20 f/z in the 
region of high temperatures. We used adimensional quantities in the figures. The temperature ranges from 10 S K to 
10 10 i^ (0.01 to 1.0 MeV). 

As an extension of this work, we are currently studying a chain of Dirac oscillators in a thermal bath. Also, 
construction of higher dimensional Dirac oscillator in finite temperature, despite its technical difficulties, would be 
very interesting and could shed light on relativistic effects in statistical and solid state physics. 

Some models used the nonrelativistic harmonic oscillator potential for describing confinement of quarks in mesons 
and baryons [1J| . More recently, some authors suggested |a, ll3( that the Dirac oscillator could be a good candidate to 
be used as the confinement potential in heavy quark systems. On the other hand, recently ultrarelativistic heavy ion 
experiments have search for the quark gluon plasma, a novel phase of QCD in which quarks and gluons are deconfined 
[l5| . We expect that studies on quark-gluon plasma models and their thermal properties could be subsidized by our 
results. 
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FIG. 3: The heat capacity for the harmonic and the Dirac oscillator for higher temperatures and TiLu/mc 2 — 1. Dashed line 
corresponds to the Dirac oscillator and the solid line the harmonic oscillator. 



ACKNOWLEDGMENTS 



M. H. P., R. R. L. and C. A. S. A. were supported in part by Conselho Nacional de Desenvolvimento Cientffico e 
Tecnologico-CNPq. 



REFERENCES 



M. Moshinsky and A. Szczepaniak, Jour. Phys. 1989 A22 L817. 
D. Ito, K. Mori and E. Carriere, Nuovo Cim. 1967 51A 1119. 
P. A. Cook, Lett. Nuovo Cimento 1971 1 419. 
H. Ui and G. Takeda, Prog. Theor. Phys. 72, 266 (1984). 

M. Moshinsky and Y. Smirnov, The harmonic oscillator in modern physics Hardwood Academic Publishers, The Nether- 
lands 1996, pp. 289-404. 

M. Moreno and A. Zentella, Jour. Phys. 1989 A22 L821. A. L. Salas-Brito, Phys. Rev. Lett. 1990 64 1643. 

J. Benitez, R. P. Martinez y Romero, H. N. Nunez- Ye pez and A. L. Salas-Brito, Phys. Rev. Lett. 64, 1643 (1990). 

P. Rozmej and R. Arvieu, J. Phys. A 32, 5367 (1999) arXiv:quant-ph/9903073 . 

Y.Nogami and F. M. Toyama, Can. J. Phys. 1996 74 114. 

F.M. Toyama, Y. Nogami and F. A. B. Coutinho, Jour. Phys. 1997 A30 2585. 
T. Loucks, Phys. Rev. 1965 A139 1333. 

R. Renan, M. H. Pacheco and C. A. Almeida, J. Phys. A 33, L509 (2000) arXiv:cond-mat/9911356 . 
F. Dommguez-Adame and M. A. Gonzalez, Europhys. Lett. 1990 13(3) 193. 
F. Ravndal, Phys. Lett. B 113, 57 (1982). 

P. Braun-Munzinger, Nucl. Phys. A 681, 119 (2001) | [arXiv:nu cl-ex /000 7021"] . 



